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We study theoretically the Rabi oscillations of a resonantly driven two-level system linearly cou- 
pled to a harmonic oscillator (vibrational mode) with frequency, Uq. We show that for weak coupling, 
ojp <C Wo, where ui p is the polaronic shift, Rabi oscillations are strongly modified in the vicinity of 
the Rabi-vibronic resonance Qr — ojo, where Qr is the Rabi frequency. The width of the resonance 
is (Qr — loq) ~ lUp^' < u>q' 13 3> uj p . Within this domain of Qr the actual frequency of the Rabi oscil- 
lations exhibits a bistable behavior as a function of Qr. Importantly, within the resonant domain, 
the oscillator is highly excited, which allows to treat it classically. 
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I. INTRODUCTION 

Coupling to environment tends to damp the Rabi 
oscillations 1 of a resonantly driven two-level system. 
Usually, the environment is viewed as medium with con- 
tinuous spectrum of modes. Less common is the situa- 
tion when environment possesses a single or several well- 
defined frequencies. For concreteness we will consider 
the situation depicted in Fig. [I] when the lower level of 
the two-level system is coupled to an oscillator (a mass, 
M, and a spring), which represents a single vibrational 
mode. Obviously, coupling to the oscillator has a strong 
effect on the Rabi oscillations in the regime of the vac- 
uum Rabi splitting^ when the oscillator frequency, ujq, is 
close to the transition frequency, u>xz. It is less obvious 
what effect the coupling to the oscillator will have on the 
Rabi oscillations when uiq is much smaller than W12 and 
is comparable to the Rabi frequency, Qr. One can argue 
on physical grounds that the effect of coupling on the 
Rabi oscillations will be strong in the vicinity of the con- 
dition, ujq R3 Qr, which we dub Rabi-vibronic resonance. 
Indeed, consider the Hamiltonian 

H = fiXhx, (1) 

describing the linear coupling. Here X — 7^75 +b) is the 

operator of the oscillator displacement, W is a creation 
operator of the vibrational quantum, hi is the occupa- 
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tion of the level E\, and fi — (2Mu)q) A, where A is 
a dimcnsionlcss coupling constant. In the course of the 
Rabi oscillations the average fi\ changes with time as 

n 1 (t) = -(l + cosQ R t). (2) 

Then at Qr rs ujq, the second term in Eq. |2]) gives 
rise to a resonant driving force acting on the oscillator. 
In turn, the stro ngly driven oscillator provides a reso- 
nant feedbadP^ on the two-level system. Thus, as Qr, 




FIG. 1: A schematic illustration of the system under consid- 
eration. Two-level system is driven by near-resonant light, 
ui\2 ~ (E2 — Ei). The level Ei is linearly coupled to a clas- 
sical oscillator with frequency ojo- The Rabi oscillations are 
strongly modified when ujo is close to Qr, where Qr is the 
Rabi frequency. 



which is proportional to the ac field driving the two-level 
system, increases, we expect the Rabi oscillations to be 
strongly modified near the resonant condition. 

Among possible experimental realizations of the situ- 
ation Fig. [l] is a suspended carbon nanotube in inhomo- 
geneous electric field, which creates a confinement for an 
excitort 5 -"® . Localized exciton can be viewed as a two-level 
system. Bending modes having discrete frequencies due 
to finite nanotube length, can be viewed as oscillators 
with very low friction. While a typical transition fre- 
quency in such a system i d 9 l 10 l w 12 ~ 10 15 Hz, the oscilla- 
tor frequency is much smaller ujq ~ 4. 3 TO 9 Hz. Resonant 
condition can be achieved by adjusting the illumination 
intensity. 

Another area in which the situation Fig.[l]is relevant, is 
the cavity QEEP3], where a two-level system is realized in 
the form of a superconducting qubit, while the oscillator 
is a LC-circuit. Majority of experimental and theoret- 
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ical studies in this field are focused on the strong cou- 
pling in the domain uq ss u±2- However, in experiments 
Refs. |12)13j an ac driven superconducting qubit was cou- 
pled to a "slow" LC-oscillator tuned to £Ir. It was ob- 
served that the noise spectrum of the oscillator exhibits 
a Lorcntzian peatP as a function of D,r — ^o- In the- 
oretical papers^H^ initiated by the experiment Ref. [12] 
collective motion of the oscillator coupled to a qubit was 
studied within the density matrix formalism, and both 
subsystems were treated quantum-mechanically. In view 
of complexity of this description, final results were ob- 
tained numerically for particular values of a coupling pa- 
rameter, A. A notable finding of Refs. [T5HT7] is that, in 
the vicinity of the condition ujq w Qr, collective Rabi- 
vibronic motion becomes bistable. 

There are still several basic questions to be answered, 
among which: 

(i) how the frequency, s, of the collective oscillations 
depends on A? 

(ii) what is the width of the resonance, i.e., the do- 
main So = Qr — ujq of the Rabi frequencies where Rabi 
oscillations are modified due to coupling? 

( in ) how the decay of the Rabi oscillations depends on 
the oscillator friction? 

The above questions are studied in the present paper. 
Our main finding is that the width of the Rabi-vibronic 
resonance is small for weak coupling, namely, 



II. BASIC EQUATIONS 

We first assume that the displacement, X(t), is a clas- 
sical variable and will later check this assumption. The 
equation of motion for X(t) reads 



X + + ulX = — (1 - w) 
' 2M y > 



(4) 



where w = 1 — 2ni is the population inversion, and 7 
is the friction in the oscillator. The equation for w(t) 
and nondiagonal elements of the density matrix have the 
form of the optical Bloch equations 



w(t) = -£l R v-T(l + w), 
«(*) = -(A-AiX(t))»-£t 

v(t) = (a - fiX(t)^u + n RW - 



(5) 



(6) 



(7) 



So = A 4 / 3 



U 



,2/3,, 1/3 



^0 < w o, 



(3) 



where T is the relaxation rate of the excited state. Note 
that, while the oscillator is driven by w(t), it exercises a 
feedback on the two-level system via u{i) and v(t). 

We require that the level E\ at t = is occupied while 
the level E2 is empty, i.e., w(0) = — 1. We also assume 
that the dipole moment and dipole current are initially 
zero, leading to v(0) = and u(0) — 0, respectively. 
From Eq. ([5| , we see that the initial conditions for v and 
w require that w(0) = 0. 



where lo p = \ 2 lo q is the polaronic shift. Eq. ^ suggests 
that, while Sq is much smaller than u>o, it is much big- 
ger than Up. Most importantly, Eq. (|3|) guarantees that, 
in the resonant domain (il^ — uq) ^do, the oscillator is 
highly excited and can be treated as classical. This allows 
the analytical description of the resonance. In this re- 
gard, the situation we consider, two-level system coupled 
to a classical oscillator, is similar to the Rabi resonance 
considered in Refs. |18)19j . where two- level system was 
driven by two classical fields: one with frequency close to 
u\2 and one with frequency close to Qr. 

We will see that in the domain (FIr — uo) ~ Sq, the fre- 
quency s of collective oscillations differs from VIr also by 
~ Sq. Bistable behavior of the dependence s(Qr) emerges 
naturally within our approach; the frequency jumps be- 
tween two stable regimes are also ~ Sq. In addition, in 
the present paper we study how the Rabi-vibronic reso- 
nance depends on detuning, A, of the driving frequency 
from U\2 , on intrinsic anharmonicity of the oscillator, and 
how the modified Rabi oscillations decay with time due 
to relaxation of the two-level system and due to friction 
in the oscillator. 



III. MODIFIED RABI OSCILLATIONS 

A. Oscillation frequency 

The system Eqs. ^-^h can be reduced to two coupled 
equations by excluding v[t) and expressing u(t) in terms 
of w(t). Then one gets 



w - 



-Tw+ln 



-Q R {A-fxX(t))u(t), 
(8) 



dt' 



«(*) = -/ ^e r ^'-^ 2 UX(t')-A) w(t')+Th+w) 

JO "R L 

(9) 

We start from the simplest case, r —> 0, A — > 0, 7 —> 0, 
and search for a solution of the system Eqs. Q, ([8]), and 
([9]), in the form w(t) — — cos st. Substituting this form 
into Eq. Q we find the displacement 



X(t) = 



/i cos st 



2M(u 



(10) 
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the ratio of the energy of oscillations to the vibrational 
quantum, u>o, we get the following estimate 



MX 2 uj$ 
co 



. 1/3 

Wo ^ A- 2 / 3 » 1. 



(14) 



Thus, for weak coupling, the classical treatment of the 
oscillator is justified. 



FIG. 2: Red line: Dimensionless frequency, z, denned by Eq. 
(17 1, of Eq. ( |19[ l versus the dimensionless deviation, x, from 
the fijj = too- Unstable solution is shown with dashed line. 
Blue line, corresponds to the absence of coupling to the oscil- 
lator. 



Static displacement, /i/2Afoj 2 , can be neglected com- 
pared to the oscillating part. Substituting X(t) into Eq. 
f9]), we find u(t) 



u{t) 



/j? (1 — cos 2st) 
'8Mn R (uj 2 - s 2 ) 



(11) 



Substituting Eqs. (10), (11) into the right-hand side of 
Eq. ((IJ), and equating the terms oc cos st in both sides, 
we find a closed equation for s 



s 2 )' 



(12) 



Thus, coupling to the oscillator causes the shift of the 
oscillation frequency from tt Rl as stated in the Introduc- 
tion. Note that the term oc cos 2st in u(t) will also give 
rise to nonresonant contribution oc cos3si in w(t), caus- 
ing a weak anharmonicity of the oscillations. Away from 
resonance, we can substitute s = tt R into the right-hand 
side of Eq. (12). Then Eq. (12) yields a correction to 



the Rabi frequency due to coupling to the oscillator 



B. Vicinity of the resonance 



To incorporate finite detuning, A, into Eq. ( 12 ) it is 



convenient to rewrite Eq. ([8j) keeping all A-dependent 
terms in the right-hand side 

w + fl 2 R w + fi 2 X(t) [ dt'X(t')w(t') - A/j,X(t) + A 2 



(15) 



A(nX(t)-A)w(t) + An / dt'X(t')w(t'). 

Jo 



The term oc A 2 in the right-hand side leads to a standard 

1/2 

modification of the Rabi frequency to (Jl R + A 2 ) . The 
last term is proportional to sin 2 st, and does not contain 
the first harmonics. The term oc cos st comes from the 
combination AfiX(t) in the left-hand side. Emergence of 
this term, which is odd in detuning, is the result of the 
coupling of the vibronic mode only to the level E\ . This 
term results in the following modification of Eq. ( 12 ) 



A — s = 



u p lJJ 



cjpCj 2 A 



«K 2 



s 2 y 



(16) 



Near the resonance (£Ir — ujo) -C ujq this equation can be 
simplified. Upon introducing dimensionless variables, 



UJQ 



2/3 1/3' 



(17) 



'0 



UJpUJQ 



(13) 



x = x 1 + 8A' ; 



Qr. — wo 

2/3 1/3 : 



(18) 



This expression is valid only if the correction on the 
right-hand side is much smaller than (coq— fijj) . Equating 
the correction to (cjq — fl R ), we find that the width of 
the resonance, (ojq — n R ) ~8q, is given by Eq. (13). 

Recall now our basic assumption that the oscillator is 
classical. We are now in position to verify this assump- 
tion. In the resonant domain the amplitude, X(t), can 
be estimated from Eq. (10) as X ~ (i/MuiqSq- Then for 



Eq. ( 16 1 assumes the form 



(z - x)z 2 + A'z 



64' 



where dimensionless detuning, A', is defined as 
A'= A 



A 1/3 2/3" 



(19) 



(20) 
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FIG. 3: Dimensionless frequency, z, of oscillations of driven 
two-level system is plotted from Eq. ( 19 1 versus the dimen- 
sionless deviation, x' , from the resonance fo r th ree positive 
dimensionless detunings, A', denned by Eq. (20 1. As detun- 
ing increases, the unstable branch shifts from positive z to 
negative z, and both stable values of z become positive (for 
positive x') or negative (for negative x'). 
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FIG. 4: The same as in Fig. [3] for three negative detunings. 
Note that negative A' broadens the range of bistability. 
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Note that characteristic detuning A ~ lu p uj 
bigger than the width of the resonance, So, but much 
smaller than f2#. Solution of Eq. ( 19 ) for zero detuning 
is plotted in Fig. [2] Blue line z = x corresponds to the 
Rabi oscillations without coupling. We see that bista- 
bility develops for x > 3 ■ 2~ 8 / 3 . At x = 3 ■ 2~ 8 / 3 , the 
frequency s experiences a jump by 2~ 8 / 3 5 - Two values of 



z corresponding to stable solutions define, via Eq. (17 1, 
two frequencies of the modified Rabi oscillations. They 
also define two corresponding amplitudes of the oscillator 

(2Afw )7 1/2 



X 



(21) 



The last factor in Eq. ( |21[ ) is the amplitude of a zero- 
point motion of the oscillator. As the dimensionless de- 
viation, x, from resonance increases, the upper branch 
approaches z = x. For this branch the frequency of the 
Rabi oscillations is close to Q,r and the amplitude of os- 
cillator is small. For the lower branch z is small, i.e., the 
frequency of the oscillations approaches luq with increas- 
ing x. For this branch the oscillator is highly excited. 

Figs. 3 and 4 illustrate the effect of detuning on the 
frequency of oscillations, s. Note that there is a quali- 
tative difference between Fig. 2 for zero detuning, and 
Figs. 3 and 4 for positive and negative detunings, respec- 
tively. For zero detuning, the domain of bistability exists 
only when O > u; , whereas for finite detuning, bistable 
regions emerge both to the left and the right from the res- 
onance. As one changes the dimensionless deviation, x' , 
from the resonance, from negative to positive, for A' = 0, 
bistability corresponds to x' > 3 ■ 2~ 8 / 3 . For finite pos- 
itive detuning, A' > 0, the first domain of bistability 
occurs at x' < 0, then disappears, and reemerges at pos- 
itive x' greater than 3 ■ 2~ 8 / 3 . Conversely, finite negative 
detuning simply broadens the domain of bistability as 
compared to A = 0. Bistable region starts for x' < 0. 
Peculiar dependence of s on the deviation from resonance 
is also reflected in the amplitude of the oscillator. This 
effect is discussed in the Sect. V. 



C. Effect of intrinsic anharmonicity of the oscillator 

Suppose that in addition to the harmonic part of 
the oscillator energy, MoJ^X 2 /2, a weak intrinsic anhar- 
monicity, kX 4 /4, is present. Then Eq. (10) will assume 
the form 



2M' 



(22) 



The second relation, s 2 = Vt 2 R — /i 2 X 2 /8, between the 
amplitude, X a , and the frequency, s, which follows from 
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Eq. ([8| remains unchanged. It is now more convenient 
to express s from this relation and substitute it into Eq. 
(22). This yields a cubic equation for X a 



2w (w - ttn)X a 



6k 



X, 



2M' 



(23) 



If we now set k = 0, then Eq. ( 23 ) will have multiple real 
solutions for X a in the domain (Tin — Ug) < —3 • 2 — 8 / 3 <5o , 
i.e., the same as determined from Eq. (19) with A' = 0. 
We see from Eq. (23) that, depending on the sign of 
k, intrinsic anharmonicity can either shift the thresh- 
old of bistability to the left (for positive k) or to the 
right (for negative k). Anharmonicity will also affect the 
magnitude of the jump of the frequency, s, of the oscil- 
lations. This magnitude will get modified from 2~ 8 / 3 <5o 

to 2- 8 / 3 <5 (l + 6K/[i 2 y 2/3 , i.e., the jump will become 
smaller for k < 0. 



IV. DECAY OF THE OSCILLATIONS 

Up to now we disregarded both mechanisms of dissi- 
pation: finite relaxation, T, and the friction in the os- 
cillator, 7. Rabi oscillations will decay with the rate T, 
which is determined by T, in the regime r > 7, or by 
friction in the regime T <C 7. We will consider both 
cases separately. We emphasize that, as the oscillations 
decay, so does the coupling between the oscillator and 
two-level system. Thus the decay will be accompanied 
by the change of frequency back to fl R . We cannot cap- 
ture this evolution of frequency with time analytically. 
To find the decay rate, T, only, we will adopt the ap- 
proach based on the energy conservation. 



A. Friction-dominated regime, 7 > T 

Upon neglecting T in Eq. ^ and setting A = we 
have 



w + n 2 R w = -fJ, 2 X(t) I dt'X{t')w{t'). (24) 
Jo 

Multiplying both sides by w and integrating from to t, 
we arrive to the following conservation law 



w 2 Q 2 R , 2 . 11 2 * r* 



db'X{t')w{t')\. (25) 



The right-hand side describes the energy exchange be- 
tween the two-level system and the oscillator. 



As a next step we multiply the equation of motion of 
the oscillator Eq. by X and integrate from to t. 
Then we arrive at the second conservation law 



X*_ 
2 



, , 2 Y 2 



+ 7 / dt'X 2 (t') 



dt'w(t')X(t'). 



At long time, t ^> T 1 we have, w 
X -> 0. Then the left-hand side of Eq 



2M 

(26) 

w — > and also X, 
251) for w turns 



to — Qfj/2. Combining Eqs. (25) and (26), we arrive to 
the relation 



n R = 2M 7 / dt'X 2 (t' 



(27) 



This relation is convenient to find the decay rate, T, be- 
cause it contains only X 2 , which is insensitive to the 
change of frequency in the course of decay. Substituting 
X 2 oc exp[— 2Yt], we find the following relation between 
T and the amplitude of the oscillations at time < 1/T 



T = - ( (MX 2 n R ) 



(28) 



Note that the second factor in the right-hand side can be 
rewritten as (MX 2 VL 2 R )/Vt R . The oscillator is classical 
when this ratio is big. Thus we conclude that T » 7, 
i.e., the Rabi oscillations in the region of resonance decay 
faster than undriven oscillator. As a next step, we distin- 
guish two cases of weak and strong friction. In the first 
case, to find X that should be substituted into Eq. ( 28 ) 
one can use Eq. ( 10 1 obtained without friction. Then 
one gets 



T = 



(uj - s) 2 



(29) 



where s is determined from the cubic equation Eq (13). 
In the region of resonance, the difference oj — s is ~~3o, 
which yields 



T 



CJ 



1/3 



(30) 



Weak friction requires that \uj ~ s\ ~ S Q » T, i.e., 7 <C 

Up. 

In the region of strong friction the difference uq — s 
should be replaced by T. Then Eq. (29) contains T in 



both sides. Upon solving this equation, we get 

jL 
u„ 



(31) 



for 7 > w p . Equations (30), (31) match when 7 ~ u p . 
The validity of this expression is limited from above 



G 



by the condition that the oscillator is classical. As 
we replace luq — s by T, the estimate for X is X <~ 
/i/(MwoT). Then the kinetic energy can be estimated 

as Mwj^/OMiT)] 2 - Wp /3 Wg /3 /7 2/3 - The condi- 
tion that it is bigger than luq limits 7 in Eq. (31) 
to 7 <C ^(wo/^p) 1 ^ 6 ! and correspondingly limits T to 

T«<5 (cV^ 1/6 - 



From Eqs. (30 1, (31) we see that, upon increasing fric- 



tion, the decay rate, T, first grows linearly with 7, and 
then sublinearly, as 7 1 / 3 . At the boundary of applica- 
bility of the classical description we have T = 7. For 
even bigger 7 classical treatment of the oscillator is not 
justified, but we expect that the oscillator will eventually 
decouple from the two-level system, and Rabi oscillations 
will proceed as they do in the absence of the oscillator. 



where S = ujq — Qr. We can now estimate m as 
flfjt c w u>otc, where t c is the time after which w^lt) be- 
comes ~ 1. Consider first the limit 8 — > 0. Then w^(t) 
grows with time as (ojpOJot/^f) . This yields m ~ 7/0J0A 2 . 
Condition that m ^> 1 should be consistent with the con- 
dition that the oscillator is classical. The latter condition 
reads: uj p ujq 3> 7 2 . In the domain when both conditions 
are met we have 1 <C m <C \- 



Consider now the limit 7 — !• 0. For St <C 1, Eq. (32) 



yields w^\t) ~ {Lo p Lo t 2 
t r . ~ (uinUJci) 12 and m ~ A -1 



This leads to the estimate 
WpWoJ and m ~ A Small-i expansion of Eq. 
( [32] ) is valid if the product St c is small. With t c found 
above, this product can rewritten in the form 



B. Initial stage of the oscillations 

After the ac driving field is switched on, the popula- 
tion inversion starts to oscillate with frequency Or. Af- 
ter some number, m, of the oscillation cycles the fre- 
quency crosses over to s. The question of interest is how 
m depends on the coupling strength, A. We will estimate 
m using the fact that at initial stage the system Eqs. 
Q, (24) can be solved perturbatively in small parame- 
ter ujp/ujQ. To find the perturbative solution, we substi- 
tute the "bare" Rabi oscillations — — cos(Qfli) into 
Eq. Q and find X(t) with initial conditions X(0) = 0, 
X{0) = 0. The obtained X(t) together with w^jt) is 
then substituted into the right-hand side of Eq. (24). 



Solving this second-order differential equation with a 
given right-hand side, we find that the amplitude of os- 
cillations becomes 1 — w^(t), i.e., it acquires a time- 
dependent correction with w^\t) given by the following 
expression 



w 



i) 



(*) 



( 7 2 
2 



7 



-AS 2 ) 2 
AS 2 



7* 
2 



(7 2 + ^ 2 ) 



7 2 + 4£ 2 



1 



-7*/2 



cos St 



(32) 



<5t r 



/J 2 \ 1 - 



(33) 



On the other hand, the oscillator can be treated as clas- 
sical when the ratio in the right-hand side of Eq. ( 33 ) is 
small. Thus, taking the limit St -C 1 in Eq. (32) is jus- 
tified, and the frequency of the Rabi oscillations crosses 
over from to s after ~ A -1 cycles. Correspondingly, 
after ~ A -1 cycles, the oscillator will "forget" about ini- 
tial phase, imposed by the initial conditions, and will 
execute a forced harmonic motion with frequency, s. 



C. Relaxation-dominated regime, f>7 



At finite relaxation rate of the two-level system Eq. 
( 26 ) assumes the form 



! 1 / p 2 \ r 2 3 



V fdt' X{t')w(t') f dt"e T{t "- tr V 2 X(t") [w{t") + T(l + to) 
Jo Jo 

(34) 

I 



Without coupling to the oscillator the right-hand side oscillations due to relaxation. Indeed, upon substitut- 
is zero, and Eq. (34) describes the decay of the Rabi ing w — flu sin(f2#i') exp(— Ti') and taking the limit 
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t — » oo, the last term in the left-hand side takes the 
value 3my8T, which leads to T = 3r/4. Natu- 
rally, this value of T follows directly from Eqs. 
(JTJ). Finite coupling to the oscillator would increase 
the decay rate only if at t —> 00 the integral in the 
right-hand side exceeds il R . Contribution of the second 
term in the square brackets to the integral can be esti- 
mated upon noticing that the product X(t")w(t") is a 
slow function. Assuming that X and w both decay as 
exp(— Tt") and that T T, the integral over t' reduces 
to J °° dt't' sin(2«0 exp(-2Tf) = T/s 3 a T/wg. Then 

one gets the estimate uj 2 TT / (ujq — s) 2 for this contribu- 
tion. Since ujq — s cannot be smaller than T, this con- 
tribution cannot exceed uj 2 which is much smaller than 
tl R . The contribution from the first term in the square 
brackets also cannot exceed Q R . This becomes appar- 
ent upon performing integration by parts after which the 
contribution from the first term assumes the form 



dt'e 



-IV 



dt"e Tt "' 2 X(t")w{t") . (35) 



If X and w decay much faster than exp(— Tt" /4), the in- 
ner integral saturates at times t' <C T. Then the contribu- 



tion Eq. (351 can be estimated as Q 2 R <->J 2 / (u)q — s) , which 



is again much smaller than VL 2 R . We thus conclude that, 
while coupling to the oscillator modifies the frequency, 
the decay of the Rabi oscillations in the relaxation- 
dominated regime is always dominated by the relaxation 
rate. 



V. NUMBER OF VIBRATIONAL QUANTA 

We studied the behavior of the of frequency of the Rabi 
oscillations near the the resonance VL R = ojq. The number 
of vibrational quanta, N, is also sensitive to the devia- 
tion, Q R — loq, from the resonance and to the detuning, 
A. This number can be expressed as follows 



N : 



1/3 



1 



16z s 



1/3 



1 



16z 2 



(36) 



where z is the solution of Eq. ( |19[ ). The dependence 
of N on dimensionless deviation, x' , and dimensionless 
detuning, A', is plotted in Fig. [5] The values of N shown 
correspond only to stable regimes of oscillations. The 
line of bifurcation points separates the (x', A') domains 
with and without bistability in Fig. [5] In the domain 
of bistability the higher and lower values of N coincide 
along the red line. Away from the red line the high-A 
and the low-A values differ very strongly. High-A values 




0.2 -1 



FIG. 5: Excitation level of the oscillator, j\f = 16iVA 2/3 , 
wher e N is the number of vibrational quanta, is plotted from 
Eq. ( 19 1 vs. dimensionless detuning and dimensionless de- 
viation from the resonance, x' = (Qr — Wo)/(k>oWp) 



The 



thick black line of bifurcations separates the "inner" domain 
of parameters (blue domain where bistability is absent) and 
outer domain (yellow domain where bistability is present). In 
the outer domain the ratio of M- values corresponding to the 
two stable solutions grows rapidly away from the boundary. 



correspond to the regime of oscillations with frequency 
close to wo, see Fig. |3j whereas low A- values correspond 
to the frequency of oscillations close to SIr. 

At x' < 0.5 and to the right from the bifurcation line 
there is no bistability. The value of A is low in this 
domain around A' = 0. As A' increases, A grows for 
both signs of detuning, A'. However, for A' > (blue 
detuning) the growth of A is monotonical. At the same 
time, for A' < (red detuning) the bistability sets in 
at certain critical A'. Upon further increase of |A'| the 
low-A value does not grow, while the high-A value grows 
rapidly. 

It is instructive to compare the results shown in Fig. [5] 
to the results of Refs. |15ll6j . The curves N(A',n R ) 
were obtained in Refs. [15116] by numerical solution of 
the system of master equations for the density matrix 
describing both the two-level system and the oscillator. 
Firstly, there is a qualitative agreement in the shape of 
the boundary of bistability. In Refs. [15116] only the 
low-A values are plotted. The prime observation made 
in Refs. [15,16] was that there is a strong difference be- 
tween these low-A values for blue and red detunings, 
namely, for blue detuning A is much higher. Our an- 
alytical results in Fig. [5] agree qualitatively with this 
observation. 
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VI. CONCLUDING REMARKS 

Frequency, flu, of the Rabi oscillations is proportional 
to the square root of the excitation power. This linearity 
has been demonstrated in many experiments. Even when 
Rabi oscillations are damped, the dependence s(£Ir) can 
be extracted from the position of maximum in the Fourier 
transform^ of the signal, w(t). We predict that, for a 
two-level system coupled to a vibrational mode, the po- 
sition of maximum of the Fourier transform will deviate 
from the linear behavior near the resonance flu = ujq. 
Both to the left and to the right from the resonance the 
position of maxim corresponds to s < Qr. The relative 
width of the resonant region depends on the coupling, A, 
to the vibrational mode as A 4 / 3 . We also predict that, 
in the vicinity of the resonance, the dependence s(CIr) 
exhibits a hysteretic behavior with two stable values of s 
corresponding to two stable regimes of the Rabi oscilla- 
tions. 

The underlying physics of the Rabi-vibronic resonance 
is the following. Without coupling, population inversion, 
w, and displacement, X, satisfy the harmonic oscillator 
equations with frequencies and uiq, respectively. With 
coupling, two-level system acts as a driving force oc w on 
the oscillator, while the back-action of the oscillator on 
w is peculiar. The structure of back-action force is ml 2 , 
as can be seen from Eqs. Q. This structure im- 

plies that back-action is of a parametrical type, i.e., X 2 
adds to fl R . Thus, at Qr ojq, it appears that Qr is 
modulated with frequency ~ 2Q,r. This, however, does 
not lead to a parametric instability. Instead, the oscilla- 
tor motion gets synchronized with the Rabi oscillations. 
In this regard, there is certain analogy to the synchro- 
nization of the Rabi oscillations to a sequence of pulses^ 
applied to the detector with repetition period chosen to 



be 2tt/Q,r. 

As it was pointed out in Introduction, the situation 
when a two-level system undergoing the Rabi oscillations 
is coupled to the oscillator is actively studied in connec- 
tion to the circuit QEEM The most common situation in 
circuit QED is when the oscillator frequency, uq, is tuned 
to the transition frequency, ui%2, of the two-level system. 
Among physical effects predicted for this domain is that 
two or multiple qubits can get strongly co upled to each 
other via coupling to a common oscillatoiP^. Rabi- 
vibronic resonance corresponds to the domain loq <C Wi2- 
Still the effects similar to those discussed in Refs. |22|23j 
will take place under the conditions of the Rabi-vibronic 
resonance. In particular, we anticipate that Rabi oscil- 
lations in two driven two-level systems with £Ir = too, 
coupled to the same oscillator will get synchronized. 

As a final remark, classical treatment of the vibrational 
mode adopted in the present paper does not allow one to 
capture the quantum jumpj^ between the stable regimes 
of collective motion of the two-level system coupled to the 
oscillator. We also did not consider the effect of thermal 
noise which leads to the activated switching^ between 
the steady regimes even within a classical description of 
the oscillator. 



Acknowledgments 

We are grateful to C. Bohme for a discussion which 
initiated this study. We thank J. Koch and F. von Op- 
pen for bringing Refs. |15ll6j to our attention. We also 
acknowledge discussions with E. G. Mishchenko, O. A. 
Starykh, and B. Spivak, and the support of the Grant 
DMR-0808842 (R.G). 



1 I. I. Rabi, Phys. Rev. 51, 652 (1937). 

2 E. T. Jaynes and F. W. Cummings, Proc. IEEE 51, 89 
(1963). 

3 L. Y. Gorelik, A. Isacsson, M. V. Voinova, B. Kasemo, 
R. I. Shekhter, and M. Jonson, Phys. Rev. Lett. 80, 4526 
(1998). 

4 Ya. M. Blanter, O. Usmani, and Yu. V. Nazarov, Phys. 
Rev. Lett. 93, 136802 (2004). 

5 C. Galland and A. Imamoglu, Phys. Rev. Lett. 101, 157404 
(2008). 

6 C. Galland, A. Hogele, H. E. Tiireci, and A. Imamoglu, 
Phys. Rev. Lett. 101, 067402 (2008). 

7 I. Wilson-Rae, C. Galland, W. Zwerger, and A. Imamoglu, 
e-print |arXiv:0911.1330| (2009) . 

8 J.-J. Li, W. He, and K.-D. Zhu, Phys. Rev. B 83, 115445 



(2011). 

9 R. B. Capaz, C. D. Spataru, S. Ismail-Beigi, and S. G. 

Louie, Phys. Rev. B 74, 121401 (2006). 
10 P. Avouris, M. Freitag, and V. Perebeinos, Nat. Photon. 
2, 341 (2008). 

1 see J. Schoelkopf and S. M. Girvin, Nature (London) 451, 
664 (2008), and references therein. 

12 E. Il'ichev, N. Oukhanski, A. Izmalkov, T. Wagner, M. 
Grajcar, H. G. Meyer, A. Y. Smirnov, A. Maassen van den 
Brink, M. H. S. Amin, and A. M. Zagoskin, Phys. Rev. 
Lett. 91, 097906 (2003). 

13 V. I. Shnyrkov, D. Born, A. A. Soroka, and W. Krech, 
Phys. Rev. B 79, 184522 (2009). 

14 A. Yu. Smirnov, Phys. Rev. B 68, 134514 (2003). 

15 J. Hauss, A. Fedorov, C. Hutter, A. Shnirman, and 



9 



G. Schon, Phys. Rev. Lett. 100, 37003 (2008). 

S. Andre, V. Brosco, M. Marthaler, A. Shnirman, and 

G. Schon, New J. Phys. 10, 095018 (2008). 

O. V. Zhirov and D. L. Shepelyansky, Phys. Rev. Lett. 

100, 014101 (2008). 

Ya. S. Greenberg, E. Il'ichev, and A. Izmalkov, Europhys. 
Lett. 72, 880 (2005). 

Ya. S. Greenberg, Phys. Rev. B 76, 104520 (2007), 
Ya. S. Greenberg, E. Il'ichev, and F. Nori, Phys. Rev. B 
80, 214423 (2009). 

D. R. McCamey, K. J. van Schooten, W. J. Baker, S.-Y. 
Lee, S.-Y. Paik, J. M. Lupton, and C. Boehme, Phys. Rev. 



Lett. 104, 017601 (2010). 

A. N. Korotkov, Phys. Rev. B 83, 041406 (2011). 

J. Majer, J. M. Chow, J. M. Gambetta, J. Koch, B. R. 

Johnson, J. A. Schreier, L. Frunzio, D. I. Schuster, A. A. 

Houck, A. Wallraff, A. Blais, M. H. Devoret, S. M. Girvin, 

R. J. Schoelkopf, Nature 449, 443 (2007). 

M. Dukalski and Ya. M. Blanter, Phys. Rev. A 82, 052330 

(2010). 

M. I. Dykman and M. A. Krivoglaz, Zh. Eksp. Teor. Fiz. 
77, 60 (1979) [Sov. Phys. JETP 50, 30 (1979)]. 



